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We consider the chemical freeze-out of H, H and (j) muhistrange hadrons within a Statistical 
Hadronization Model inspired approach. We study particle yields across a wide range of reaction 
energy and centrality from NA49 at SPS and STAR at RHIC experiments. We constrain the physical 
conditions present in the fireball source of strange hadrons, and anticipate results expected at LHC. 
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Introduction: We study multistrange hadron pro- 
duction in the context of the quark-gluon plasma (QGP) 
formation in relativistic heavy ion collisions Given 
the relatively small reaction cross sections of multistrange 
hadrons in hadron matter, the observed yields of S(gss), 
S, f2(sss), f2, (f>{ss) [l-Q are considered probes of the 
earliest stage of the QGP-fireball hadronization. 

The yields of these particles have been considered pre- 
viously within a global approach, see e.g. Q. Here we 
show that it is possible to analyze multistrange hadron 
yields alone. When this is done we find that multistrange 
and non-strange hadrons share the same freeze-out con- 
dition. We will discuss the meaning of this discovery 
below addressing dynamics of hadronization. We also 
address the forthcoming LHC effort to measure multi- 
strange hadron yields in high multiplicity pp Q , and soon 
after, in A+A reactions. 

QGP hadronic particle production yields are usually 
considered within the statistical hadronization model 
(SHM) SHM has been successful in describ- 

ing (strange) hadron production in heavy ion collisions 
for different colliding systems and energies. These re- 
sults showing successful global fits of particle yields in 
the SHM framework inspired us to study multistrange 
hadron yields alone in this separate analysis for the pur- 
pose of: i) establishing that SHM is appropriate for de- 
scribing yields of these particles, ii) assessing if their 
yields are consistent with the established bulk matter 
properties of the QGP fireball, thus testing the single 
freeze-out hypothesis for particles with large and small 
hadron reaction cross sections, and iii) understanding 
better how the future LHC results may help arrive at 
a distinction between SHM model approaches. 

SHM Models: We begin by introducing the three 
principal SHM approaches: 

a) Taking the view that SHM has a limited theoretical 
foundation, one can seek simplicity in an effort to obtain 
a qualitative description of the yields for all hadrons with 
just a small number of parameters. An additional attrac- 
tion is that this assumption leads to a model with chem- 
ical equilibrium hadron yields is explored. The main re- 
sult of this approach is that the hadronization in high en- 
ergy heavy ion collisions at RHIC requires T > 175 MeV, 
and this high value is close to the lattice crossover tem- 
perature, between deconfined and hadron phase [l3l.[T^. 

b) In order to arrive at a precise description of the 



bulk properties, such as strangeness and entropy content 
of the hadron fireball, we need precise capability to ex- 
trapolate hadron yields to unobserved kinematic domains 
and particle types. This is achieved by introducing sta- 
tistical occupancy parameters 7^ > l,i = q,s. Within 
this approach there is good systematic behavior of phys- 
ical observables as a function of collision conditions such 
as energy or centrality 8, 15 -19]. The yields of hadrons 
are in general found not to be in chemical equilibrium, 
7i ^ 1; the hadronization temperature is found near to 
T ~ 140 MeV. 

While this value of T could be further away from 
the deconfinement crossover domain, this is where chi- 
ral symmetry restoration is achieved [l^, [2^, and QGP 
is transformed into hadrons. Moreover, in this approach 
the variation of the freeze-out temperature with bary- 
ochemical potential parallels the slope seen in the lattice 
data. Another important outcome of this approach is 
that a fit to data offers a good statistical significance. 
Results obtained can be interpreted in terms of a dy- 
namical picture of nearly chemically equilibrated QGP, 
decaying into free streaming hadrons. The high intrin- 
sic QGP entropy content explains why equilibrated QGP 
turns into chemically overpopulated (over-saturated) HG 
phase space - the fast breakup of QGP means that the 
emerging hadrons do not have opportunity to re-establish 
chemical equilibrium in the HG phase. 

c) Single Freeze-out or/and Strangeness Nonequilib- 
rium model has as the main objective statistically signif- 
icant description of hadron yields achieved with minimal 
effort. Only strangeness chemical non-equilibrium is al- 
lowed. This is often enough to produce a decent data fit 
and to assure that all particles can be formed at the same 
physical condition pl| - |25| . The main result of this ap- 
proach is a hadronization temperature near T ~ 160 MeV 
which agrees with Hagedorn temperature [1^, [l^l • 

Particle ratios of interest: We must include in our 
theoretical consideration of multistrange hadron yields 
the contributing yield of decaying hadron resonances. 
Within SHM these individual yields generally depend 
on several parameters. The phase space occupancy 7^ 
scales particle yields according to light quark content, 
and a similar parameter 7s refers to strange quark con- 
tent. Temperature T quantifies the size of accessible 
phase space. Baryo-chemical potential /ub differentiates 
baryons from antibaryons and strange chemical potential 
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/is does the same for strangeness. There is also a poten- 
tial /ii3 related to different number of up and down quarks 
which is constrained by proton and neutron asymmetry 
in colliding nuclei, and the overall yield is normalized by 
a volume parameter V. 
By considering the ratio 



we eliminate in good approximation most of the SHM 
parameter dependencies since; 

i) By taking the product of particle and antiparticle, we 
eliminate baryo-chemical potential /iB as well as strange 
chemical potential fis- 

ii) We also eliminate the strange quark phase space oc- 
cupancy 7s, because the strange and anti-strange quark 
content in the numerator and denominator is the same, 
iii) The overall normalization is eliminated by the fact 
that we have the same number of hadrons in the ratio 
numerator and denominator. 

The S/(/) ratio depends on the probability of finding 
a non-strange d, d-quark at the formation of S~((iss) 

and S^((iss), respectively. This is expressed by the light 
quark phase space occupancy jq. Furthermore, temper- 
ature T controls the magnitude of 

the (non-relativistic) phase space ratio of and 4>. 
S(1321) is always a decay product of S*(1530). Thus 
aside of the ground state i = 1 : S(1321), 51 = 2 one 
must include in the sum the 5*(1530), 92 = 4 resonance. 
Consideration of this special yield ratio parallels the ear- 
lier effort made to identify ^s/lq in Ref. 2§]. 

We extend our considerations to include single 
strange K"'"(us), K~(us) mesons and triple strange 
VL~ {sss),VL [sss) baryons considering the ratios: 

Given the quark content, both S/K and ^l/cj) are pro- 
portional to strange quark yield, i.e. the strange quark 
phase space occupancy 7s and a function fi{T). 

The arguments leading to Eq. ([J) , Eq. ([3]) are strictly 
valid only in Boltzmann approximation. Considering 
quantum statistics, there is some residual dependence 
of f{T) on chemical parameters, involving higher pow- 
ers of 7g for the ratio S/c/) Eq. ([1]), and higher powers of 
7s for the ratio S/K Eq. In order to estimate the 
magnitude of the quantum statistics effect we calculate 
the actual particle ratios with SHAREv2 [ll| using both 
quantum and Boltzmann statistics. We find that the 
Boltzmann approximation we used overestimates S/0 by 
0.25%, which is always negligible. For S/K, we find that 
it is overestimated by Boltzmann approximation by up 



TABLE I: Fit parameters used to determine particle yields 
for incompatible centrality bins using /(A'^part) = a - A^pa^t +c 
(see text for details). 





a 


b 


c 


TT 


4.179 X 10"' 


1.072 


7.107 X 10"' 


7V+ 


4.247 X 10"' 


1.048 


6.422 X 10"' 


K+ 


5.433 X 10"^ 


1.111 


-1.014 X 10"' 


K- 


4.812 X 10"^ 


1.107 


-3.859 X 10-2 




1.228 X 10"^ 


1.247 


-4.678 X IQ-^ 




8.978 X 10-* 


1.221 


9.390 X 10-'' 


4>" 


4.162 X 10"^ 


1.203 


-9.311 X 10-^ 



to 4%, the relatively larger effect is due to the relatively 
low mass of the kaon. Since the experimental error is 
much greater we continue to consider the simple theo- 
retical Boltzmann yields. When further below we con- 
sider ratios involving pions, all results are obtained using 
SHAREv2, which accounts for resonance decays and all 
yields can be obtained using quantum statistics. 

Experimental data: We consider An data from the 
CERN-SPS NA49 experiment, and for the STAR ex- 
periment at RHIC the acceptance rapidity interval is 
\y\ < 0.5; therefore at RHIC we use the yield per unit 
of rapidity dN/dy and omit the differential dy when re- 
ferring to relative yields. For the (f) meson we consider 
the recently published data from STAR fl| and the up- 
dated data from NA49 [3|. We collected the necessary 
data for S and S baryons from Refs. 4-6]. 

We do not use NA49 158 GeV results, since these ex- 
perimental results do not allow to interpolate the differ- 
ent centrality bins used to measure different multistrange 
particle yields. We could not simply combine data from 
different centrality bins seen the variation of yields with 
centrality (that is 7s). The STAR 62 experiment provides 
data in several centrality bins, defined as a percentage of 
the most central collisions: data from the most central 
collisions is found in the centrality — 5% interval and the 
most peripheral collision results presented are in 70—80% 
bin. The relation to A^part and/or impact parameter b is 
discussed in 0. 

We use recent data for K* mesons from STAR exper- 
iments at both y/SNN = 200 and 62.4 GeV from [?[. For 
the SPS NA49 data we use yields from [H, H^. 

We note that different centrality bins are often chosen 
for different particle types. In order to be able to form 
particle ratios in a common centrality interval, we in- 
ter/extrapolate, that is fit individual yields as a function 
of the number of participants using a simple functional 
form /{Npart) = c • iVpg^j-t -t- c. We show the fit parameters 
a,b and c in TableUand compare the experimental results 
and the fit in figure [T] 

Particle ratios: After this preparation we can form 
ratios of particle yields as shown in figure [2] and table 
im We note that the S/(/) relative yield does not change 
much over a wide range of energies and centralities, in 
contrast to the individual hadron yields which enter the 
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TABLE II: Values of ratios 2/0 Eq. ([T]), H/K Eq. ([S]), H/tt and (jy/Ti Eq. Q obtained from the data and the resulting estimated 
uncertainty in 7^ and 7^, respectively. When symbol 'E' is shown in the error column, the data ratio is result of interpolation 
and/or extrapolation needed to account for different centrality bins. 



Experiment 


Centrality 


X 10 


5lq 


H/K X 10^ 




H/tt X 10'^ 


(t>/K X 10 


(t>/Tl X 10^ 


STAR 62 


0-5% 


3.04 


E 


4.19 


9.6% 


6.22 


1.38 


2.04 


QT'AT? cry 


c; 1 no/ 


o.UU 


111 


A HQ 
4.U0 


y.z/o 


D.zU 


l.oD 


Z.UD 


STAR 62 


10-20% 


2.94 


E 


4.06 


9.3% 


5.98 


1.38 


2.04 


STAR 62 


20-40% 


2.88 


12.5% 


3.79 


E 


5.48 


1.32 


1.91 


STAR 62 


40-60% 


2.85 


14.6% 


3.38 


E 


4.65 


1.18 


1.63 


STAR 62 


60-80% 


2.49 


19.3% 


2.84 


E 


3.45 


1.14 


1.38 


STAR 200 


0-20% 


3.02 


11.8% 


4.06 


12.9% 


6.04 


1.34 




SPS 80 AGeV 


7% 


3.33 


24.5% 


3.04 


22.7% 


2.60 


0.83 


0.88 


SPS 40 AGeV 


7% 


2.45 


42.1% 


1.89 


18.0% 


3.23 


0.78 


0.83 


SPS 30 AGeV 


7% 


2.57 


66.5% 


1.85 


24.3% 


2.10 


0.63 


0.72 


" For STAR 200 (j) 


/tt considering 


figure 14 in 


2] we give an 


average of data for centralities up to 50%. 
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FIG. 1: (Color online) Data points (full symbols) of particle 
yields used in the analysis, and their respective fitted central- 
ity dependence. 



FIG. 2: (Color online) Data points of H/</> Eq. ©, H/K 
Eq. (O, H/tt and (j)/-K Eq. (g}. The straight line for H/0 = 
0.281 . 



ratio. The average value of all available data points is 
S/(/) = 0.281 with an error at 15% level. 

The remarkable result, the constancy of E!/(/) means 
that at SPS and RHIC energies the mechanisms and con- 
ditions at which double-strange particles are produced 
are very similar and that, according to Eq. ([1]), there is 
a constraint between values of 7g and T, which we now 
explore in figure |3] where we show in the T, 7^ plane the 
theoretical SHM results as lines for a constant ratio S/i/). 
These values are obtained using SHAREv2 and varying 
7^ and T, with all other model parameters fixed to a rea- 
sonable physical values. In this way we also confirm once 
again the analytical formula Eq. ^ . 

We limit the magnitude of 7^ by a critical value of 
light quark phase space occupancy 7^''''. 7^0 = 7^ < 
(7q"*)^ = exp(mjro/T), which is the condition where 
the pion phase space distribution function diverges for 
m^rO = 135 MeV/c^. The experimental values S/0 ~ 
0.281±15% are found consistent with all SHM models in 
that for 7, = 1 we find the value T = 170 ± 10 MeV, and 



for 7, ^ 1.63 a value T ^ 140 MeV. 

Test of SHM models: We have seen that hadroniza- 
tion of S and (j) is consistent with the three different SHM 
models, but there is an interesting constraint between 
T, 7g arising from the constancy of the relative S and (j) 
yield. We see also in figure [2] that the variation of S/K 
is significant, it changes by a factor of 2.3 in the entire 
data range. Considering that we already have established 
by the study of S/(/) that the hadronization temperature 
does not vary this indicates that there is a variation of 
7s value by a factor of about 2.3 in the data range. We 
conclude that a fixed value 7s = 1 cannot be chosen. 
This rules out the SHM model a) . We also note that this 
argument can be made in the same way considering the 
variation of the other ratios in fig. [5J e.g. S/tt and (p/ K. 

SHM results for S/K and fi/t/f in T, 7s plane are shown 
in figure |4l obtained by the same method as before; i.e. 
using SHARE with other SHM parameters fixed at an 
appropriate value. For a given S/K and/or $7/0 a slight 
7g dependence remains, since there are unrelated reso- 
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120 130 140 150 160 170 180 190 

T [MeV] 

FIG. 3: (Color online) Lines of a constant given ratio Eq. IT]) 
H/0 e [0.03,0.45] in the r,7q plane. The lines for 7, = 1 
and 7q = 7,"^'' are presented by solid black lines. The average 
result, 0.281, of all SPS and RHIC experiments is highlighted 
by a thick gray line. As this ratio is considered constant, this 
line indicates the prediction of LHC results. 

nances decaying into K (and to lesser degree S). Thus 
we present for each fixed value of S/if two extremes: 
7g = 1, and 7^ — 7"'*. The effect is depicted in figure |4^ 
in terms of two lines shown by the same line type. 

Note that similarly as for 7g, there is a critical value 
for 7s based on the Bose-Einstein condensation of the rj 
meson (r/ = 0.55(uu + dd) + 0.45ss [H 113). The large 
values of 7s could be relevant to the future LHC results. 
To compare theory and experiment we show the thick 0- 
20% STAR 62 line and by looking at the bottom frame 
of figurelHwe obtain the prediction: 5.5 x 10~^ < 51/0 < 
7.0 X 10~^, the variation due to variability of hadroniza- 
tion temperature. 

In order to further elaborate the validity of models b), 
c) we show on the right in table|lT]the ratios S/tt, 0/K, 
and 0/7r, where 



7r~V7i-7r+' K~Vk-K+' 7r~V7r-7r+' ^' 

The experimental S/tt and i^/tt relative yields vary by 
a factor ~ 3.5 in both cases. In figure [5] we show the 
(jj/iT ratio and compare to theory as a function of 7s at 
fixed given T. Model b) with T ~ 140, 7, = 7crit im- 
plies that the different experimental results correspond 
to 1 < 7s < 2.4. These values are consistent with the 
large value of 7^ = 7crit — 1-6. On the other hand, for 
7q = 1 several fixed T lines nearly coincide in the inter- 
esting range 210 > T > 160 MeV. This means that the 
growth in yield of </> is nearly compensated by the growth 
in TT-multiplicity. It will be very interesting to see how 
LHC results will line up in this presentation, since we see 
that the high energy RHIC results even at 7^ = 1 imply 
7s > 1. A value 7s > 1 is incompatible with the picture 
of strangeness production in hadron collisions, and im- 




120 130 140 150 160 170 180 190 
T [MeV] 

FIG. 4: (Color online) Lines of constant ratio H/K (a) and 
Q,/4> (b). Experimental data from most central 0-20% STAR 
62 are indicated by a thick line (b), and is assumed in the 
bottom frame as a prediction. See text for more detail. 

plies the presence of a strangeness dense QGP phase as 
a source of hadrons. 

Behavior at LHC: As already remarked the ratio 
S/0 ~ 0.28 is firmly constrained and cannot change. 
Even under the extreme LHC conditions we expect that 
this ratio will be the same as at RHIC. However, consid- 
erable changes can be expected for other (multi) strange 
particle ratios which were discussed earlier [13, Ha] . Here 
we will mainly address the ^/tt ratio. 

When we accept the premise, that entropy and 
strangeness are conserved during the hadronization, we 
can predict values of the phase space occupancy 7s in 
chemical semi- and non-equilibrium models for LHC. We 
expect a 20% increased value of strangeness over entropy 
s/S ~ 0.037 [11]. For the two models under considera- 
tion (T = 140 MeV/^ = ^critical g^^^ MeY,-/q = 1) 
this value suggests [SJ] "/s/"/q — 1-55. The expected (jj/ir 
ratio is 2.95 x 10~^ and 3.90 x 10~^ for the two models as 
indicated by the boundaries of the LHC band in figure [S] 
Experimental results of this magnitude requires 7s > 1 
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FIG. 5: (Color online) The relative (ji/vr Eq. @ yield as a 
function of 73 in several hadronization scenarios, see text. 
The vertical solid black line shows the chemical equilibrium 
with 7s = 1. For experimental data see table |IT] Predicted 
values for LHC are indicated in blue. 

and concludes in favor of chemical non-equilibrium the 
still ongoing discussion of chemical equilibrium models. 

Summary and conclusions: We find that the rel- 
ative particle yield S/(/) is practically constant as func- 
tion of centrality and reaction energy at RHIC and SPS. 
We find that these particles, despite their small reac- 
tion cross-sections are emerging at the same hadroniza- 
tion condition as all bulk particles. This result was an- 
ticipated 33] for a fast expanding QGP fireball which 
under-cools and rapidly breaks apart (hadronizes), and 
has been used extensively in single hadronization mod- 



els [H [H [H HI] 

Variation in the ratio S/K (and thus also (/>/K cx 
7s /7g) implies a variation in strange phase space oc- 
cupancy 7s, in agreement with the expectation that 
strangeness production grows with energy and centrality 
of the collision. This experimental result is incompatible 
with the chemical equilibrium model a), for which also 
the parameter 7^ is fixed to 1 by definition. 

Considering further the yields S/tt and ^/vr consis- 
tency with the bulk matter particle production rates is 
arrived at within the chemical non-equilibrium model b) 
with 7q > 1 and 7s > 1. These values imply that the 
observed strange hadrons yields are above chemical equi- 
librium, a feature predicted to be signature for hadroniza- 
tion of a QGP-fireball [S^. The expected further increase 
of 7s > 1 at LHC implies a further increase of the 4>/tt 
ratio, providing a clear distinction between chemical non- 
equilibrium model b) and semi-equilibrium model c). 

Our results show that the yields of all multistrange 
hadrons available today arc 1) compatible with the SHM 
picture of hadron formation, 2) are well described by cur- 
rent chemical nonequilibrium hadronization models in 
the parameter domain obtained from the other hadron 
yields, 3) these data are incompatible with the chemical 
equilibrium single-freeze-out SHM. A critical test of our 
approach is that in LHC-ion experiments the S/^ ratio 
remains the same as has been observed at SPS and RHIC. 
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